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^ ! Abstract 

The measurement of forward jet cross sections has been suggested as a promis- 
ing probe of new small x parton dynamics and the question is whether the new 



CO \ HERA data provide an indication of this. In this paper the influence of re- 



solved photon processes has been investigated and it has been studied to what 
extent the inclusion of such processes in addition to normal deep inelastic scat- 
. tering leads to agreement with data. It is shown that two DGLAP evolution 

Q^-i! chains from the hard scattering process towards the proton and the photon re- 

Qh| spectively, are sufficient to describe effects observed in the HERA data, which 

have been attributed to BFKL dynamics. 

> 
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1 Introduction 

Experimental data from deep inelastic scattering (DIS) in a kinematic region where new 
parton dynamics is expected to become noticeable, i.e. at small values of the scaled proton 
momentum, x, are not described by models based on interactions with pointlike photons. 
In a previous paper |lj we have demonstrated that the addition of interactions through 
resolved photons offers a possible explanation of the observed discrepancies and leads to 
good agreement with all available data. 

This paper is devoted to a more detailed discussion of the resolved photon concept and 
comparisons with data on forward jet production in DIS, since the forward jet cross section 
has been advocated as a particularly sensitive measure of small x parton dynamics 0, 
Analytic calculations based on the BFKL equation 0, || in the leading logarithmic 
approximation (LLA) are in fair agreement with data. However, recent calculations of 
the BFKL kernel in the next-to-leading logarithmic approximation (NLLA) || have given 
surprisingly large corrections, and it remains to be shown whether the data can still be 
reasonably described. 

Monte Carlo generators based on direct, point-like photon interactions (DIR model), 
calculated from leading order (order a s ) QCD matrix elements, and leading log parton 
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showers based on the DGLAP evolution do not take any new parton dynamics in the 
small x region into account and are therefore not expected to fit the experimental data. 
Recent results from the HI J7| and ZEUS pi experiments on forward jet production 
exhibit significant deviations from the predictions of such models. Also, next-to-leading 
order calculations (NLO, i.e. order a 2 ) assuming point-like photons predict too small a 
cross section compared to data. 

The study of forward jet production with contributions from direct as well as resolved 



photon processes has been performed using the RAPGAP 2.06 ]9|, [L0| Monte Carlo event 
generator. 



2 Resolved Photons in DIS 

In electron-proton scattering the internal structure of the proton as well as of the ex- 
changed photon can be resolved provided the scale of the hard subprocess is larger than 
the inverse radius of the proton, 1/-R 2 ~ Aq CD , and the photon, 1/R 2 ~ Q 2 , respectively. 
Resolved photon processes play an important role in photo-production of high jets, 
where Q 2 ~ 0, but they can also give considerable contributions to DIS processes (JTTJ, O 
if the scale fi 2 of the hard subprocess is larger than Q 2 , the inverse size of the photon. 
This led to the idea of including contributions from resolved photon processes as part of 
deep inelastic scattering to get a more complete description fl"3|j . 

In the following we give a brief description of the model for resolved virtual photons 
used in the Monte Carlo generator RAPGAP. Given the fractional momentum transfer 
of the incoming electron to the exchanged photon, the Equivalent Photon Approximation 
provides the flux of virtual transversely polarized photons ]10|, |9], and references therein]. 
The contribution from longitudinally polarized photons has been neglected. The partonic 
structure of the virtual photon is defined by parameterizations of the parton densities, 
x 7 / 7 (x 7 , /i 2 , Q 2 ), which depend on the two scales ji 2 and Q 2 [TJ, |TB . The following 
hard subprocesses are considered (RES model): gg — > qq, gg — > gg, qg —> qg, qq — > gg, 
qq — > qq, qq —>■ qq. Parton showers on both the proton and the photon side are included. 
The generic diagram for the process q^g p — > qg including parton showers is shown in 
Fig- 0- 

Since the photon structure function depends on the scale, /i 2 , of the hard scattering 
process, the cross section of resolved photon processes will consequently also depend on 
the choice of this scale. It has to be carefully considered in which range of /x 2 ' IQ\ the 
photon-parton cross section can be factorised into a parton-parton cross section convoluted 
with the parton density of the photon. The parton density of the photon is evolved 
from a starting scale Q\ to the scale fi 2 , the virtuality at the hard subprocess, giving a 
resummation to all orders. 



2.1 Parton Distribution Functions 

Due to factorization of the cross section, the parton densities of both the virtual photon 
and the proton enter into the calculations. The proton structure function, F 2 , has been 
measured to high accuracy and therefore the various parameterizations only give marginal 
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Figure 1: Deep inelastic scattering with a resolved virtual photon and the q^g p — > qg 
partonic subprocess. 



differences in the measurable kinematic region. Two parameterizations of the parton 
distribution in the proton, GRV 94 HO (DIS) and CTEQ4D, have been considered, which 
both give good agreement with the proton structure function data |[7], |18|| . It was found 
that the results produced were identical at the percent level, when keeping Aqcd fixed. 
In the following we use only GRV HO. 

The photon can interact via its partons either in a bound vector meson state or as 
decoupled partons if the pt of the partons is high enough. The splitting 7 — > qq is called 
the anomalous component of the photon. The structure function of virtual photons has 
been measured only recently, but not nearly to the same precision as the proton structure 
function. However, it turns out that data are in good agreement with the parameterization 
of Schuler and Sjostrand (SaS) fll5f . The SaS parameterization offers a choice of Ql values 
at which the anomalous part becomes effective. We have studied these choices resulting 
in different magnitudes of the parton densities, and consequently of the cross sections. 
For the SaS parameterization we have used Ql as given by eq.(12) of ref. |y|] (IP2 = 2). 
This choice is also suitable for a description of other hadronic final state properties (not 
considered in this paper), like energy flow, forward particle spectra and jet cross sections. 

The hadronic contribution to the virtual photon structure function decreases rapidly 
with increasing Q 2 , which means that the main contribution at large Q 2 comes from the 
anomalous piece in the photon splitting. This is completely calculable in pQCD and leads 



to an expected agreement between the parameterization of Gliick - Reya - Stratman [14 
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and that of Schuler - Sjostrand [15], but also with the simple ansatz of Drees - Godbole 



T6fl . However differences exist in the way the hadronic part of the structure function is 



matched to the pointlike part, which just reflects the theoretical uncertainty. 



In this study we will restrict ourselves to the SaS parton distributions fL5 |. 



2.2 Choice of Scale 

In leading order a s processes, the renormalization scale /j,r and factorization scale fiF 
are not well defined, which allows a number of reasonable choices. There are essentially 
two competing effects: a large scale suppresses a s (/i 2 ) but gives, on the other hand, an 
increased parton density, x/(x,/i 2 ), for a fixed small x value. The net effect depends on 
the details of the interaction and on the parton density parameterization. 

In previous papers |19|, [l] we have tried different scales like /i 2 = and /i 2 = Q 2 +p T , 



and found that these choices gave similar results. 

However, in resolved virtual photon processes, the choice of the scale /i 2 , at which the 
photon is probed, is severely restricted pt| . In a partonic process a + b — > c + d, where 
a, b, c, d denote four-vectors and where parton a has the virtuality Q 2 , the transverse 
momentum p\ of parton c is given in the small angle limit (— t <C s) by: p\ = s(—t) /(§ + 
Q 2 ), with s and t being the usual Mandelstam variables. In a t channel process the 
virtuality is given by ji 2 = —t. Thus we have:Q 



fi i = -i = p 2 T + Q 2 -^<Q 2 +p 2 T (1) 



From eq.(0) we see that the scale fi 2 is always larger than the transverse momentum 
squared of the hard partons and less than Q 2 +p\. In the following we will use fi 2 = Q 2 +p\ 
as scale for both resolved virtual photon processes and for direct photon processes. This 
choice of scale provides a smooth transition from the kinematic region of normal DIS into 
the range where resolved photons start contributing and further into the photo-production 
region. The same scale has also been used in NLO calculations including resolved photons 



in deep inelastic scattering pTl E2 



A basic test that the scale is reasonable is that the parton shower evolution scheme 
should not be able to produce partons with transverse momenta larger than those pro- 
duced by the matrix element for the hard scattering process. 



3 Forward Jets 

HERA has extended the available x region down to values below 10 -4 , where new parton 
dynamics might show up. Based on calculations in the LLA of the BFKL kernel, the cross 
section for DIS events at low x and large Q 2 with a high p\ jet in the proton direction (a 
forward jet) || |||] is expected to rise more rapidly with decreasing x than expected from 
DGLAP based calculations. New results from the HI and ZEUS || experiments have 
recently been published. The data can be described neither by conventional DIR Monte 

1 We are grateful to T. Sjostrand for pointing out this simple explanation 
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Carlo models nor by a NLO calculation, while comparisons to analytic calculations of the 
LLA BFKL mechanism has proven reasonable agreement. 

It should be kept in mind that both the NLO calculations and the BFKL based 
calculations are performed on the parton level whereas the data are at the level of hadrons. 
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Figure 2: The forward jet cross section as a function of x for pt jet > 3.5 GeV 
(a) and pr jet > 5 GeV (6) as measured by HI [|7| . Also shown are the RAPGAP 
predictions for the sum of direct and resolved processes (solid line) as well as the 
resolved photon contribution alone (dashed line). The RAPGAP predictions for the 
sum of direct and resolved processes without initial and final state parton showers 
(solid line), including only final state (dashed line) and only initial state parton 
showers (dotted line) are shown in (c). 
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In Fig. and b the forward jet cross section as measured by the HI collaboration p] 
is compared to the prediction of the RAPGAP Monte Carlo generator for both resolved 
photon process alone (labeled RES) and for the sum of direct and resolved processes 
(labeled DIR+RES). The calculation is performed with a scale /i 2 = Q 2 + p\ used in 
the parton densities and in the determination of a s for both direct and resolved photon 
processes. Only if both processes are included, is the Monte Carlo able to give a good 
description of the measurements. The forward jet cross section as measured by the ZEUS 
experiment || can be equally well described by this Monte Carlo program using the same 
structure functions and parameter setting. 

Recently a full NLO calculation of the forward jet cross section was performed [p^j 
with contributions from both direct and resolved virtual photons. Comparisons with HI 
and ZEUS data exhibit good agreement while it is found that LO calculations fall below 
the experimental data. The assumption is that the parton showers in the RAPGAP model 
account for higher order effects, which would explain the good description of data from 
this model. 

As a consequence, the influence of the initial and final state QCD cascade on the 
forward jet cross section has been studied in more detail and compared to the calculations. 
In Fig. we show the contribution to the forward jet cross section of HI |/J coming 
from pure matrix element calculations and the effect of including initial and final state 
radiation. It is observed that the final state radiation gives essentially no additional 
contribution to the matrix element cross section whereas the initial state radiation plays 
an important role in bringing the Monte Carlo predictions to into agreement with data. 
We, thus, come to the same conclusion as in P^fl , that the LO order DIR+RES matrix 
elements are not enough to describe the measured forward jet cross section. However, with 
initial and final state parton showers included, the data are described well and it seems 
that the parton showers are able to simulate contributions from higher order processes. 
We also find that the LO cross section calculations of [22| agree well with the RAPGAP 
results in LO order, i.e. when parton showers are switched off. It should be noted that the 
contribution from the direct process with initial and final state parton showers included, 
corresponding to the difference between the RES+DIR and RES histograms in Fig. |2]c, is 
too small to give agreement with data. 

It is also important to make sure that the hardness of the scale doesn't lead to parton 
cascades which produce radiation with larger p T than produced by the hard scattering 
matrix element. In Fig. [3] the ratio of the transverse momentum of any initial state parton 
q\ to the transverse momentum p\ of the hard scattering process is shown in the 7*p CMS 
for events which satisfy the forward jet analysis criteria. The solid line corresponds to 
ji 2 = Q 2 + the dotted line to /i 2 = 4 • p\ and the dashed line to fi 2 = p\. One can see 
that essentially all partons coming from the initial state cascade have transverse momenta 
smaller than the partons of the hard scattering p?, which is expected in a DGLAP type 
evolution, where the transverse momenta are ordered in qx towards the hard scattering 
process. Thus we conclude that the scale /x 2 = Q 2 + p T which has been used for this study 
fullfils the requirements of a DGLAP type initial state cascade. 

In Fig. ^ the different contributions of the total resolved photon cross section are shown 
separately within the cuts of the forward jet analysis. From Fig. it is observed that 
the hadronic part of the virtual photon structure function, as expected, gives a negligible 
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Figure 3: The ratio R = q\jp\ °f the transverse momenta q\ of partons from the 
initial state cascade to the transverse momentum p\ of the partons from the hard 
scattering process. The solid line corresponds to the scale fi 2 = Q 2 +p?, the dotted 
line to fi 2 = 4 • and the dashed line to fi 2 = p\. The distribution is normalized 
to the total number of entries. Please note the logarithmic scale on the y - axis. 

contribution to the measured cross section, since it dies off rapidly with increasing Q 2 . 
Fig. fp? shows that the subprocess q^g p — > qg contributes the most to the resolved photon 
cross section in the forward jet region (~ 60%) and that the subprocesses qq — > qq, 
qq — > qq and gg — > gg each give a contribution of the order of 10%. 

A small fraction of the DIS events, fulfilling the selection criteria for forward jets, 
actually contains two identified jets. Analytic calculations (in LLA) |23| have been per- 
formed in the same kinematic region and with the same jet selection as defined for the HI 
one-jet sample. The predicted ratio varies from 3% to 6% as x increases from 0.5 ■ 10~ 3 
to x = 3 • 10~ 3 . Our previously reported prediction from the RAPGAP generator 
including both direct and resolved photon processes was that about 1 % of the total for- 
ward jet sample contains two forward jets. This is about a factor of 3 lower than the 
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Figure 4: Different contributions to the total cross section of resolved virtual photons 
within the cuts of the forward jet analysis 0. In (a) is shown the ratio R = - ^ t t , 
i.e. the anomalous (pointlike) part (solid line) and the hadronic part (dashed line), 
respectively, of the resolved virtual photon cross section divided by the the total 
resolved photon cross section as a function of x. In (b) is shown the ratio R = — — — 

v ' Ores, tot 

as a function of x for different subprocesses i: qg — > qg (dotted line), qq — > qq and 
qq — > qq (upper solid line), gg — > gg (dashed line), qq — > gg (dashed-dotted line) 
and gg — > qq (lower solid line). 

prediction from the BFKL calculations but a large part of this discrepancy could be due 
to hadronization effects which would reduce the prediction of the parton level BFKL cal- 
culation. Recently this ratio has been measured by the HI experiment to be 1 %, in 
excellent agreement with the prediction of RAPGAP. This gives further confidence in the 
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basic concept of resolved photons even at large Q 2 . 

The ZEUS collaboration has presented a measurement of the forward jet cross section 
as a function of E^/Q 2 in the kinematic region Q 2 > 10 GeV 2 , y > 0.1, E e > > 10 
GeV, r] jet < 2.6, x jet > 0.036, E T jet > 5 GeV, p zbreit > GeV, 2.5 ■ 10" 4 < x < 8 • 10~ 2 
but without implementing the DGLAP suppression cut 0.5 < E\jQ 2 < 2. The results 
were compared to the predictions from different Monte Carlo programs and the conclusion 
is that only the RAPGAP DIS generator including interactions through resolved virtual 
photons can describe the data over the full range in E^ /Q 2 ■ In Fig. the ZEUS results are 
shown together with the prediction of RAPGAP and in Fig. |5|6 the contributions coming 
from the hadronic and the pointlike part of the virtual photon structure function are 
presented separately. The influence of initial and final state radiation is shown in Fig. 
which again illustrates that the initial state radiation gives an essential contribution to 
the cross section while the contribution from final state radiation is negligible. 



4 Summary and Discussion 

Recent experimental data on forward jet production show deviations from traditional LO 
Monte Carlo models assuming directly interacting point-like photons. It is tempting to 
assume that the observed effects could be explained by BFKL dynamics. 

In the present study we have shown that the addition of resolved photon processes to 
the direct interactions in DIS leads to good agreement with the data. This agreement 
does not depend on any specific choice of scale or tuning of any other parameters in the 
RAPGAP generator. The best evidence of the universality of this approach is that, with 
the same parameter setting, it is possible to describe a wide range of other data like 
the transverse energy flow [25], transverse momentum spectra of single particles ||26|| , the 



(2+1) jet rate |27j and single inclusive jet cross sections ||11|| , as we have shown in 

We have observed that the dominant contributions to the resolved photon processes 
come from order a 2 diagrams with the hard subprocess q 1 g p — > qg (see Fig. [I]). Since the 
partons which form the photon remnant per definition have smaller px than the partons 
involved in the hard scattering, a situation with non qx ordering is created. 

In the LO DIR model, the ladder of gluon emissions is governed by DGLAP dynamics 
giving a strong ordering of q? for emissions between the photon and the proton vertex. The 
models describing resolved photon processes and BFKL dynamics are similar in the sense 
that both lead to a breaking of this ordering in q?. The BFKL picture, however, allows 
for complete dis-ordering in qj>, while in the resolved photon case the DGLAP ladder is 
split into two shorter ladders, one from the hard subsystem to the proton vertex and one 
to the photon vertex, each of them ordered in q t (see Fig. [l]). Only if the ladders are long 
enough to produce additional hard radiation might it be possible to separate resolved 
photon processes from processes governed by BFKL dynamics. Thus the resolved photon 
approach may be a "sufficiently good" approximation to an exact BFKL calculation and 
the two approaches may prove indistinguishable within the range of x accessible at HERA. 

It should be emphasized again that the usual NLO calculation assuming point-like 
virtual photons contains a significant part of what is attributed to the resolved structure 
of the virtual photon in the LO scheme. PTJ. The NLO calculations including contribu- 
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tions from direct as well as resolved photons, however, are similar to using LO matrix 
elements for both the photon and the proton with the addition of parton showers, as im- 
plemented in the RAPGAP generator. Recent calculations have proven good agreement 
with the predictions of RAPGAP and thus also with data. This indicates that higher 
order contributions are well simulated by the inclusion of parton showers in LO Monte 
Carlo generators. 
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Figure 5: The cross section of forward jets as a function of E^/Q 2 as measured 
by ZEUS [p4fl . In (a) the data are compared to the prediction of RAPGAP. The 
solid line shows the sum of direct and resolved virtual photon contributions, whereas 
the dashed line shows the resolved photon contribution alone. In (b) is shown the 
part of the cross section coming from the anomalous component (solid line) and the 
one coming from the hadronic component of the virtual photon separately. In (c) 
the RAPGAP predictions are shown for the sum of direct and resolved processes 
without initial and final state parton showers (solid line), including only final state 
(dashed line) and only initial state parton showers (dotted line). 
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